The gravitational energy-momentum flux 



g ■ J. W. Maluf * and F. F. Faria 

(N ; Institute de Fisica, 

^! Universidade de Brasilia 

< : CP. 04385 

\0 '■ 70.919-970 Brasilia DF, Brazil 

; and 

^ ! K. H. Castello-Branco 

^ ■ Institute de Fisica, Universidade de Sao Paulo 

O : Caixa Postal 66318, 05315-970 Sao Paulo SP, Brazil 

, 

o : 

m , 

. Abstract 
O . 

We present a continuity equation for the gravitational energy- 
momentum, which is obtained in the framework of the teleparallel 
• • . equivalent of general relativity. From this equation it follows a gen- 

I eral definition for the gravitational energy-momentum flux. This def- 

^ ■ inition is investigated in the context of plane waves and of cylindrical 

. Einstein- Rosen waves. We obtain the well known value for the energy 

flux of plane gravitational waves, and conclude that the latter exhibit 
features similar to plane electromagnetic waves. 
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1 Introduction 



Attempts to define the gravitational energy in tetrad theories of gravity were 
first put forward by M0ller [1], who noticed that the tetrad description of the 
gravitational field allows a more satisfactory treatment of the gravitational 
energy-momentum. He observed that a suitable expression constructed out 
of the tetrad field and of the torsion tensor could possibly yield a definition 
for the gravitational energy density. The torsion tensor cannot be made to 
vanish at a point by a coordinate transformation. This fact refutes the usual 
argument against the nonlocalizability of the gravitational energy, and which 
rests on the reduction of the metric tensor to the Minkowski metric tensor 
at a point in space-time by means of a coordinate transformation. 

It is well known that the tcleparallel equivalent of general relativity 
(TEGR) [2, 3, 4, 5, 6, 7, 8, 9, 10, 11] provides an alternative description 
of Einstein's general relativity, in which the gravitational field is described 
by the tetrad field. In fact the very first attempt to construct a theory of the 
gravitational field in terms of a set of four linearly independent vector fields 
in the Wcitzcnbock geometry [12] is due to Einstein [13]. 

A well posed and mathematically consistent expression for the gravita- 
tional energy has been investigated [14] . It arises in the realm of the Hamilto- 
nian formulation of the TEGR [15] and satisfies several crucial requirements 
for any acceptable definition of gravitational energy. The definition of the 
gravitational energy-momentum is obtained from the constraint equations 
of the theory. The results so far obtained indicate that such expression is 
physically relevant. 

In any small neighborhood of space the gravitational field can be consid- 
ered constant and uniform. The principle of equivalence asserts that in such 
neighborhood it is always possible to choose a reference frame in which the 
gravitational effects are not observed (we are adopting Einstein's version of 
the principle of equivalence [16]). Thus in such reference frame we should not 
detect any form of gravitational energy. Therefore it is reasonable to expect 
that the localizability of the gravitational energy depends on the reference 
frame, but not on the coordinate system. In fact any other form of relativistic 
energy depends on the reference frame. It turns out that the gravitational 
energy definition presented in Ref. [14] displays the feature discussed above, 
namely, it depends on the reference frame. More precisely, it depends on 
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the choice of a global set of tetrad fields since the energy expression is not 
invariant under local SO (3,1) transformations of the tetrad field, but is in- 
variant under coordinate transformations of the three-dimensional spacelike 
hypersurface (reference frames are better conceived in terms of fields of vector 

bases [17]). 

It is pointed out in Ref. [14] the existence of a preferred global frame that 
allows a proper discussion of the energy, momentum and angular momentum 
of the gravitational field, and that will be briefly recalled in Sec. II of this 
article. In this section we will also present an alternative interpretation of the 
tetrad field, based on the relation between inertial and noninertial frames. In 
Sec. Ill we consider the Lagrangian field equations of the TEGR and obtain 
the continuity equation for the gravitational energy-momentum. From this 
equation it follows a definition for the gravitational energy-momentum fiux. 
The latter is applied in Sees. IV and V to plane waves and Einstein-Rosen 
waves, respectively. By means of simple calculations we obtain the well 
known value of the gravitational energy fiux carried by plane gravitational 
waves, and we show that this value coincides with the fiux of gravitational 
momentum along the propagation direction, in similarity with plane electro- 
magnetic waves. 

Notation: space-time indices fi, u, ... and S0(3,l) indices a, b, ... run from 
to 3. Time and space indices are indicated according to /i = 0, i, a = (0), (i). 
The tetrad field e"^ yields the definition of the torsion tensor: T"^,^ = 
d^e"- u — d^e"' ^j,. The flat, Minkowski space-time metric is flxed by rfab — 
ea^efe,^'^'^^ (- + ++). 



2 An interpretation of the tetrad field 

For a given space-time metric tensor g^j^^ the tetrad field is defined by the 
relation 

e"^e^i.?7a6 ^ Qnv , (1) 

where r]ab is the fiat Minkowski metric tensor. For a metric tensor g^i, there 
exists an infinit set of tetrad fields that satisfy the relation above, and which 
are related by a local Lorentz transformation, ef" ^{x) = A"- b{x)e^ ^{x). At 
every space-time point the S0(3,l) matrices A^^ satisfy 
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h.\A^dr]ab = Vcd 



(2) 



Tetrad fields are necessary in order to establish the coupling of Dirac spinor 
fields with the gravitational field [18, 19]. 

Let us consider the construction of tetrad fields for the fiat Minkowski 
space-time. For the latter we have the relation 

^e^ ^r]ab = rji^iy ■ (3) 

We note that Eq. (3) for the tetrad field is similar to Eq. (2) for the 
matrices A" Therefore the tetrad field in fiat space-time can be considered 
as a transformation from an inertial frame with coordinates x'^ to another 
inertial frame with coordinates such that dq"^ — e" ^dx^ and e"' ^ — d^q"-. 
In this case it is formally possible to integrate dq"" = e° ^rfx'^ in order to 
obtain = q"'{x^) everywhere. Because of this property the trasformation 
x'^ — > q"- is called holonomic. 

In the TEGR [14] every manifestation of the gravitational field other than 
the fiat Minkowski space-time is characterized by a set of tetrad fields such 
that df^e"- 1, — d^e"- ^ ^ 0, in which case the tetrad field is no longer given by 
gradient functions of the type ej^ = df/f. In this case the transformation 
dq"' = i^dx^ cannot be globally integrated, and the transformation is called 
anholonomic. Therefore a nontrivial manifestation of the gravitational field 
is characterized by a nonvanishing torsion tensor, T" = 9^e" ,^ — dj^e"- ^y^O 
[14]. The four-dimensional space with coordinate differentials dq"', endowed 
with Minkowski's metric tensor, will be denoted the reference space-time. 

Let the tetrad field be such that T" ^ Q. In any sufficiently small 
neighborhood of the space-time the transformation Ax'* ~ Ca^^q", where 
Ca ^ are the inverses of e" ^, may be interpreted as a transformation from a lo- 
cal inertial frame with coordinate differentials Ag" to a noninertial one with 
curvilinear coordinate differentials Aa;^. Therefore in such small neighbor- 
hood of the space-time the tetrad field measures the deviation of the physical 
space-time from a hypothetical fiat space-time. From this point of view, the 
principle of equivalence is naturally built into the kinematical description 
of the space-time, since it establishes the local equivalence of accelerated, 
noninertial frames with the gravitational field. 

We consider again the flat Minkowski space-time. In Sec. IV of Ref. [14] 
it is discussed the six conditions on the tetrad field such that the reference 
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space-time with coordinates q'^ is neither related by a boost transformation 
nor rotating with respect to the physical space-time with coordinates x^. 
These conditions uniquely fix the tetrad field. In cartesian coordinates they 
are given by 

e« ° = . (5) 

Let us consider a transformation between two cartesian space-time coor- 
dinates and q"" = q"'{x^) such that q"" is rotating with respect to x^. The 
tetrad field e" ^ = 9^^" will acquire off diagonal components in the spatial 
sector {e(j)j} [14]. Therefore the imposition of Eq. (4) prevents the rotation 
between the two cartesian space-times. Another transformation of general 
character between and is a Lorentz boost. It is easy to see that the 
imposition of the time gauge condition e(j) ° = = e*^°) k prevents the two 
space-times from beeing related by a boost transformation [14] . 

Conditions (4) and (5) will be assumed to hold also for an arbitrary space- 
time metric tensor. They have proved to be essential in the evaluation of the 
irreducible mass of the Kerr black hole [14]. Therefore for a given metric 
tensor g^^ they estabhsh a preferred global frame that allows a suitable dis- 
cussion of the energy, momentum and angular momentum of the gravitational 
field. 

3 The continuity equation and the gravita- 
tional energy-momentum flux 

The vacuum Lagrangian density of the TEGR is given by 

= -A;eS«^X,,, (6) 
where k — l/(167r), E"'"^ is defined by 
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Y^abc \^ahc _|_ rj-}hac rpcab^^ _|_ vf'^T^^ (7) 

and T"" — T'^ ij The quadratic combination Tf^'^Tabc is proportional to the 
scalar curvature -R(e), except for a total divergence [8]. 
The vacuum field equations read 



It is possible to prove by explicit calculations that 



]-e{Ra^{e)-]-ea^R{e)\ , (9) 



^ga/. 2 I "'"^ ' 2 

The quantities on the right hand side of Eq. (9) are constructed out of the 
curvature tensor i?a6/ii/(e). 

The gravitational energy-momentum is obtained in the framework of the 
Hamiltonian formulation of the TEGR [15]. The constraint equations of 
the theory are interpreted as equations that define the energy, momentum 
and angular momentum of the gravitational field. The gravitational energy- 
momentum is defined by [14] 

P"^ - I d^xdjU."^ , (10) 
where 11"* is the momentum canonically conjugated to Cai and reads 



+ e"^(^°^r'= - /^T° ij) - 2{gV - g"'e^'^)T^ . (11) 

Equation (10) satisfies all requirements that are expected from any definition 
of gravitational energy. This definition: (i) vanishes for the Minkowski space- 
time, in which case the tetrad field is given by gradient functions of the type 
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d^q"", (ii) yields the ADM (Arnowitt-Deser-Misner) and Bondi energy 



in tiie appropriate limits [20]; (ill) yields the appropriate value for weak and 
spherically symmetric gravitational fields; (iv) yields the irreducible mass of 
the Kerr black hole [14]; (v) is invariant under coordinate transformations 
on the three-dimensional spacelike hypersurface, and transforms as a vector 
under the global S0(3,l) group. 

We proceed now to derive the continuity equation for the gravitational 
energy-momentum. The field equations (8) can be rewritten by multiplying 
it by the inverse tetrad fields. We obtain 



Restricting A to spatial components, i.e., X — j, we find 



(12) 



-9o(-4A;eE«o^) - 9,(-4A;eE«*^^) = -kee'^^iAE^'^^n,^ - 5lT!"^n,a) ■ (13) 

Note that S"^'^ = — S"'^^. Taking the divergence of the equation above with 
respect to j yields 

-9,-9o(-4A;eE"'^^) = -A;a,[ee"'^(4E''^^r;,c^ - 5iT!"^''n,a)] , (14) 

or 



By integrating Eq. (15) in a space volume V we arrive at 



(15) 



d 
di 



- / d^xd.W^ 



V 



-k [ rf5,[ee"^(4S^'^'r,,^ - ^^E^^'^T,,,)] . (16) 



The time derivative of the gravitational energy in a volume V of the three- 
dimensional spacelike hypersurface is define to be minus the gravitational 
energy-momentum flux 



d 
di 



d^x diW"^ 



V 



(17) 



where 
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s 



(18) 



Wc will give one step further and assume that Eq. (18) is defined also 
for open surfaces S. For later purposes we define the gravitational energy- 
momentum fiux density (j)"'^ , 



which represents the a component of the flux density in the j direction. 

In Sees. IV and V we will carry out concrete applications of Eqs. (18) and 
(19). Here we just mention that for tetrad fields with appropriate boundary 
conditions, in an asymptotically fiat space-time, Eq. (16) leads to the total 
conservation of the gravitational energy-momentum. 

Before closing this section we briefly mention that Eq. (19) is similar to 
the gravitational gauge current defined in Ref. [21], but the interpretations 
of (p^^ are different in the two situations. In the latter Ref. the quantity 
similar to Eq. (19) is taken to represent the energy and momentum of the 
gravitational field, whereas in the present context it really represents the flux 
of energy and momentum. 



4 The energy-momentum flux of linear plane 
gravitational waves 



We will consider initially the simplest realization of a linear plane wave, 
that is a solution of Einstein's equations. The energy carried away by these 
waves is known in the hterature [22]. As in the latter Ref., we will restrict 
considerations to just one polarization of the wave, which is described by the 
metric tensor 



ds^ = -dt^ + [1 + /+(t - z)]dx^ + [1 - f+{t - z)]dy^ + dz^ . (20) 



where (/+)^ << 1. The tetrad field that satisfies conditions (4) and (5) and 
that yields the metric tensor above reads (lines are labelled by a and columns 



(19) 



by /^) 
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e"" n{t,x,y,z) 



fl 0\ 

VTTU 



Vo 1/ 



(21) 



The determinant e = det(e"^) is given by e = ^/T'^^J^^/Y^^J^. The con- 
travariant components of the metric tensor (which is given in diagonal form) 
are g'' = -1, g'' = 1/(1 + /+), g'' = 1/(1 - /+), g^' = 1. 

The nonvanishing components of the torsion tensor are easily calculated, 



f 



(1)01 



Ti 



(2)02 
T(i)i3 
2^(2)23 — 



/ 

/; 



where the dot and the prime denote derivatives with respect to the co- 
ordinates t and z, respectively. Prom the components above we evaluate 



T: 



e'^xTat.u- We find 



1 



101 



202 



13 



223 



2/+ ' 



(22) 
(23) 
(24) 
(25) 



We first calculate the energy flux . It is given by 
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s 

= -4A; / dSj e^°°e(°) o(S^^^rioi + 2^2^7202) ■ (26) 

In order to evaluate the tensor T,!^"^ we need the traces T'^ = g^'^T^ — 
g^^T^x^. We find 

T^^T^^O, = ^''^''Tiia + g^^g^^T223 (27) 

By considering j = 1 in Eq. (26) we obtain 



E221r7-i rp 
J- 202 — J- 202 



1 1 

_^221 _|_ 21221 _ 21122 j _|_ _^^21ji2 _ ^222^1^ 



0, 



because T^^i = ^22^22^11^1^^^ ^ Similarly, for j = 2 in Eq. (26) we 

have 



v112t^ 

Zj J 1 



101 — -^101 



1^^1112 _|_ 21II2 _ ^1211^ _|_ Ij^^lSjil _ ^ll2^2-j 

4 2 



. 



Therefore the only contribution to comes from the integration on 2; = 
constant surfaces, namely, on the surface orthogonal to the propagation of 
the wave. We have 



$(0) = -4A;|^d53e^°V'^)o(S^''Tioi + E2^^r202) 
= Ak j^dS^e^] 



1 1 
_ (^2^113 _|_ jiiis _ jisiij _|_ _j-^i32ii _ ^117^3^ 



1 j^2i223 _|_ 2i223 _ 2i322^ _|_ \^^2Zrjp2 _ ^22^13 j 



+ 7^202 

= -2k I 
Js 

By substituting Eqs. (22-25) into Eq. (28) we find 



dS?, e 



g g g (2101^^223 + 7202rii3) 



(28) 



^-k [ dSs , ■ (29) 

Since the function /+ is such that << 1, the equation above can be 
rewritten to a very good approximation as 

/ dS^Uf^. (30) 

s 

Assuming now that the function /+ can be written in terms of an amph- 
tude A and a frequency u; as [22] 

f+{t- z) = Acosuj{z-t) , (31) 

we have 

/+ = Au sina;(;2 ~ t) ^ f+~ sina;(;2 — t) , 

from what follows 

^^'^^^^in'u;{z-t) [ dS, , (32) 
167r Js 

where S is an arbitrary z = constant surface, and we have substituted k = 

l/(167r). We take the avarage value < ^^^^ > over a period T by considering 

the integral 



L 



T J" 

dt sin^(a;i — cuz) — — . 



Taking into account the integral above in Eq. (32) we obtain the average 
flux density < 0*^°^^ > per unit period T, flowing along the z direction, 

f " ~327 ■ ^'^'^> 
This is precisely the value obtained in the literature by means of a completely 
different analysis of the the energy flux of plane, linearised gravitational waves 
[22]. 

By evaluating the momentum flux components we arrive at a quite 
interesting result. For this purpose we first need to calculate the expression 
of the scalar E^'^'^Tbcd, which is easily obtained as 
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S^^'^T,,, = -2{g'W'T^,sT223 + 9'WT,oiT202) • (34) 
The momentum flux is given by 

+ A; dS2[4ee^''^'^^''%^.l] + k dS3[4ee^'^'T.''''%^.i] . (35) 
After some simple calculations we obtain 



^'•"X.! = -li9''g''g''T,o,T202 + 9'VV'T^,,T22s) , (36) 



and 



w''X,,^j:^^%,,^o. (37) 

By substituting Eqs. (34), (36) and (37) into Eq. (35) we conclude that 
$(1) = 0. 

A similar result is obtained for , which reads 

$(2) = k f dS^[ee^'''^''{AW^X,2-5iY.''^nJ)] 
J s 

= k j^dS, [4ee(2)2E'^^^T^,2] + k J^dS2[ee^'^^^E''''%.2 - T,'''''n^)] 
+ k [ dSs[Aee^^^^W'%,2] ■ (38) 
In this case we have 

^'"'X^2 = -\{9''9''9''T,o,T202 + 9''9''9''T,,,T22:^ , (39) 

J:^''%,2 = ^^''%u2 = . (40) 
By considering Eqs. (39) and (40) we conclude that Eq. (38) reduces to 

$(2) = 0. 
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Finally, for ^^^^ we have 



s 

S J s 

+ k [ dSs[ee^'^\4E^-'%,s - E^'^'^n^)] . (41) 
s 

Simple calculations yield 

= ^''''%.3 = , (42) 

and 

^^''%u3 = -9''9^^9^^Tn3T223 . (43) 
Substitution of Eqs. (34), (42) and (43) into Eq. (41) leads to 

= k j dS3[ee^''^ 3g''i2g''g''g''T^oiT,o2 - 2g'' g'' g''TmT22s)] , (44) 

s 

With the help of Eqs. (22-25) we find 

$(3) = k I dSs^^^^^^lM^ , (45) 

•Is 2^1 - iuy 

where S is now an arbitrary z — constant surface. Taking into account Eq. 
(31) we observe that {f\.f + {UT = (A + f+? - V'J+ = -2/;/+, since 
/; = -/+. Thus we find 

$(^) ^-k( dSs . ^ -A; / /+/; , (46) 



in similarity to Eqs. (29) and (30). 

Therefore the expression for is identical to the one for We 
obtain 

$(3) ^ ^ sin2 ^(^ f . (47) 
12 



Taking again the avarage value < > over a period T we obtain average 
momentum flux density < ^'-'^-'^ > per unit period T, flowing along the z 
direction, 

We conclude that the energy and momentum flux density per unit time 
of a linear plane wave, along the direction of propagation, are the same. It 
follows that the fluxes satisfy the relation 

^''^^nab = . (49) 

We note that a similar relation must be satisfied by the energy-momentum 
four-vector of massless particles. We also remark that plane electromagnetic 
waves display the same feature, namely, the fluxes of energy and momentum 
are the same, in natural units (c=l). 



5 The energy flux of Einstein- Rosen waves 

Einstein-Rosen waves arise as exact solutions of Einstein's field equations 
[23]. They describe cyhndrical waves determined by two functions 7(p, t) 
and ip{p,t). In cylindrical coordinates the waves are described by 

ds' = e^(^-^) (-di^ + dp') + p'e-'^d<p' + e'^dz' . (50) 
The functions 7 and ip satisfy 

^" + 1^' _ ^ = , (51) 
P 

i = pW? + {^?] , (52) 

7 = 2pV/V', (53) 

where now the prime denotes differentiation with respect to p. 

The expressions for the gravitational energy contained within a cylindri- 
cal region of arbitrary length L and radius p, around the z axis, and the 
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corresponding energy flux are very simple, as we will see. The procedure for 
obtaining these quantities is also very simple. First we determine the tetrad 
field that satisfies conditions (4) and (5), and that leads to Eq. (50). It is 
given by 



/AO 0\ 

Acoscf) — pC sin (f) 

A sin pC cos (f) 

VO BJ 



(54) 



where 



A = e^~^ , B = e^, C = e'^ . (55) 

It follows that e = det(e"^) = pe'^^'^~^\ The nonvanishing components of 
the torsion tensor are 



T{Q)Q\ = A' 



T{i)Qi = i cos , T(i)02 = -pC sin , r(i)i2 = {A-C - pC') sin , 



T(2)oi = A sin , T(2)02 = pC cos , T(2)i2 = - ( A - C - pC) cos , 



(3)13 



which in turn lead to 



Tool = AA, Tioi = AA, T202 = p'CC, 



T212 = pC{C + pC" - A), T303 = BB, T313 = BB' . 

In order to calculate the gravitational energy associated with the metric 
tensor (50) we can either calculate the a = (0) component of Eq. (10) or, 
more simply, take advantage of an equality that holds in the time gauge 
condition [14], 



14 



- / d'^xdiU^''^' = — I d^xdiieT) . (56) 
Jv Stt Jv 

The field quantities on the left hand side of the equation above are defined 
on the space-time, and those on the right hand side are defined on a three- 
dimensional spacelike hypersurface E. We note that the tetrad field (54) 
satisfies the time gauge condition because e^°^ fe = 0. In view of Gauss theo- 
rem we define the energy contained within a surface S by [24] 

where e and are calculated on E. Therefore we consider triads e(j)j on a 
t — constant surface, 

(A cos 4> —pC sin \ 
Asin(/) pCcos(/) , (58) 
b) 

and construct Qij and the inverses g'^^ and e^*^-' on the spacelike hypersurface 
E (because the metric is diagonal, g'^^ on E and on the space-time are the 
same). In this case the determinant e = det(e(j)j) is given by e = pe'^~^, and 
the trace = g^^Ti = g^\g^^T22i + g^^Ts^i) reads 

yi^ ie-'(^-^)(e^-l) . (59) 
P 

We also have T"^ = = 0. The gravitational energy enclosed by a cylinder 
of length L and arbitrary radius p is obtained by a trivial integration of Eq. 
(57), 

E{t,p) = ^e-^^-^\e''-l). (60) 

The energy per unit length e for very small values of 7 and ip reads e{t, p) ^ 
7/4. In this limiting case the latter value coincides with the value obtained 
by Thorne [25] in the analysis of the C-energy of Einstein-Rosen waves. 

We proceed now to calculate the gravitational energy flux through the 
cyhndrical surface of length L and radius p. In view of Eq. (18) it is given 

by 
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= -4k [ dS,ee^'^\E'^^%02 + ^^'%03) 
J s 

= -4k [ dS, eeW°(E22ir202 + ^^^^^303) ■ (61) 

s 

In the evaluation of Eq. (61) wc arc taking into account definitions (54) and 
(55). After simple calculations we obtain 

E^^^r202 + E^^^Taoa = ■^9^^g^^9^^{T2i2T303 + ^3132202) 

- l{9''9''9''V,,2 + 9''9''9''no,no3). (62) 

By substituting the expressions of T^ux we eventually find 

- -2k I d(j>dz [e-(^-'^)(e^ - 1)^] . (63) 
Considering a cylindrical surface of length L and radius p we obtain 

$(0) ^ ^Le-i-r-^) (eT - 1)^; . (64) 

The simplicity of expressions (60) and (64) is an indication that the 
present framework is suitable for discussing the energy-momentum properties 
of the gravitational field. 



6 Discussion 

We have derived a simple expression for the energy-momentum flux of the 
gravitational field. This expression is obtained on the assumption that Eq. 
(10) represents the energy-momentum of the gravitational field on a volume 
V of the three-dimensional spacelike hypersurface. Of course the consistency 
of the present results may be taken as a further, a posteriori justification for 
Eq. (10) to represent the energy and momentum of the gravitational field. 
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Application of Eq. (18) to linear plane gravitational waves shows that the 
latter has properties similar to plane electromagnetic waves: the energy and 
momentum fluxes are the same, in natural units. This is the main result of 
this article. Moreover, the value of the energy flux obtained here is exactly the 
same one derived in the literature by considering the energy supplied by the 
waves to a nearly continuous distribution of oscillators on a plane orthogonal 
to the direction of propagation of the wave [22]. In the latter analysis it is 
evaluated the reduction of the amplitude of the gravitational wave as it passes 
through the conflguration of oscillators describe above. Einstein-Rosen waves 
were also analysed. The expressions of the gravitational energy contained 
within a cylinder around the z axis and the corresponding energy flux through 
the cylindrical surface turned out to be simple as well. Altogether, the results 
described above indicate the consistency of the present framework in the 
investigation of the energy-momentum properties of the gravitational field. 
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